Superconductor insulator transition in thin films driven by an orbital parallel 

magnetic field effect 
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We study theoretically orbital effects of a parallel magnetic field applied to a disordered super- 
conducting film. We find that the field reduces the phase stiffness and leads to strong quantum 
phase fluctuations driving the system into an insulating behavior. This microscopic model shows 
that the critical field decreases with the sheet resistance, in agreement with recent experimental 
results. The predictions of this model can be used to discriminate spin and orbital effects. We find 
that experiments conducted by A. Johansson et al. are more consistent with the orbital mechanism. 
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Introduction and results - Applying a magnetic field 
to a disordered superconducting (SC) film, can drive it 
into a strong insulating (INS) state. This was observed 
both when the film isplaced perpendicular to the field, in 
InO 0,0, i, MoGe Q, TiN @ and NbSe 0, or when the 
film is in parallel to the field orientation in InO 0, S] and 
Bi . While several theoretical models support a SC - 



INS transition in the perpendicular orientation [lfl lllj , 
the mechanism that drives the transition when the field 
is parallel to the film, remains unclear. 

We study analytically a microscopic model of a disor- 
dered SC film with magnetic field applied parallel to the 
film, focusing on the induced orbital effects. We show 
that the field reduces the stiffness of the SC phase, lead- 
ing to strong quantum phase fluctuations manifested as 
an INS behavior. We find that the transition to the INS 
phase occurs at a critical field that depends on the SC co- 
herence length, £o, the film's thickness, d, and the sheet 
resistance, Rq [see Eq. ((TJ]. We will show that this re- 
lation does not depend on the detailed mechanism that 
drives the transition, that it allows to determine experi- 
mentally if spin or orbital effects are dominant in parallel 
field, and that the measurements of A. Johansson et al. Q 
are more consistent with the orbital mechanism. 

The emergence of a SC-INS transition in disordered 
films induced by a perpendicular magnetic field is consis- 
tent with several theoretical scenarios. A perpendicular 
magnetic field penetrates the film in the form of vortices. 
As the field increases, these vortices were predicted to 
delocalize and Bose-condense leading to an INS behav- 
ior [n| ■ An alternative numerical work studied the effect 



of thermal phase fluctuations in a disordered 2D SC [11 1 . 
The perpendicular magnetic field was shown to destroy 
phase correlations between SC islands. Conversely, a 
parallel spin-exchange field causes the order parameter 
phase and amplitude to vanish abruptly. While existing 
theories can account for the qualitative behavior seen in 
the perpendicular field orientation, they do not explain 
the surprisingly similar observed phenomenology when 
the field is parallel to the film @, @, 0| • 



We study the previously disregarded orbital effect of a 
parallel magnetic field applied to a disordered SC film. 
We find that the field uniformly decreases the SC order 
parameter and reduces its phase stiffness. As a result, 
quantum fluctuations of the phase and amplitude are en- 
hanced and can lead to an INS behavior. Our main pre- 
diction is that the critical field, B c that marks the onset of 
the INS behavior depends on the critical temperature T c , 
the film's sheet resistance, Rq and thickness d as 



B 2 JH 2 = l/2{\n(R Q /R n )-\n2K c }. 
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is the flux quantum, vq is the density of states, Rq = 
h/(4e 2 ) is the resistance quantum, Aq(B) is the mean 
field order parameter in the presence of a pair break- 
ing field [12j, and K§ is the critical value of the stiff- 
ness coefficient. Three points should be stressed herein. 
First, while the INS behavior can be a result of cither 
the proliferation of topological phase excitations, or of 
strong gaussian phase fluctuations, the detailed mecha- 
nism that drives the transition will merely change the 
numerical factor Kq, as long as dissipation in the cores is 
negligible [3, Second, we do not consider the effects 
of a spin-exchange field, as some of the materials that ex- 
hibit the SC-INS transition, such as MoGe, NbSe and Bi, 
are expected to have strong spin-orbit scattering, which 
would smear spin polarization effects. However, while the 
spin mechanism is expected to depend on the thickness, 
d, only implicitly through T c , the orbital critical filed is 
inversely proportional to d. Hence, by studying the de- 
pendance of B c field on the parameters of the system 
such as Rq, d and T c , one can determine which mecha- 
nism dominates the transition. Finally, we find that the 
data of Ref. @ is more consistent with Eq. ([1]), see Fig. Q] 
To gain insight to how the reduced stiffness can lead 
to an INS behavior we consider as an example its ef- 
fect on quantum vortex excitations, corresponding to 
a vortex loop in the D = 2 + 1 dimensional system. 
The energy cost of a circular vortex loop is -Eioop ~ 
E COIC L/C + EjL/£lnL/£, where E COTe is the energy per 
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unit length to destroy the superconductor in the core of 
the vortex and Ej is the energy per unit length to rotate 
the SC phase, and is determined by the phase stiffness. 
The logarithmic divergence will be cutoff for more com- 
plicated loop shapes. The entropy of a loop is determined 
by counting all possible configurations. On a cubic lattice 
this can be estimated by Si oop ~ ln(2D — 1) L ^ ~ L/£. 
A parallel magnetic field will reduce the core and rota- 
tional energies of the vortex loops [see discussion follow- 
ing Eq. ([3]) and j5])], leaving the entropy unchanged. As 
a result, vortex excitation become increasingly favorable 
and proliferate at a critical field which, in the absence 
of dissipation inside the vortex cores 13 , [3] marks the 
onset of the INS behavior. In addition to its effects on 
topological excitations, the reduced phase stiffness also 
enhances gaussian fluctuation of the SC phase. Loss of 
phase rigidity due to strong gaussian fluctuations give 
similar estimates for the transition field as in Eq. (fTJ), 
differing by a numerical factor Kq. 

Fig. [1] compares B c calculated using Eq.JT]) with the 
data of Ref. |a The theoretical curve was plotted with 
Kq and vq used as fitting parameters. The microscopic 
model used to obtain Eq. ffl is valid for B c /H c ii < 1. 
The clean sample of Ref. |8[ exhibit a transition into 
a metal at high magnetic fields. We identify B c in the 
cleanest sample with H c \\ . While our theory is not ap- 
plicable for the clean sample, we use this point to de- 
termine the relation between H c n, T r and i?n- The re- 



sulting H c \\ = 1.60o^/ S^^tQ 1S larger than the critical 
field calculated in Ref. by a factor 1.6. This may be 
a result of the finite thickness of the films used in Ref. @. 
The inset of Fig. [1] shows the scaling of the measured B c 
with the critical temperature T c , as expected for a spin 
mechanism. We find that the data of Ref. H is more con- 
sistent with Eq. JTJ, see caption of Fig.[TJ We have tried 
to fit the experimental data of Ref. |9| to our model. B c 
obtained using Eq. (JT|) is in better agreement with the 
experimental data than the critical field obtained from a 
linear B c oc T c or a square root dependence B c oc \fT~c-, 
as naively expected from a spin mechanism, with weak 
or strong spin-orbit scattering, respectively. However, as 
apposed to InO [|| , the distinction between the two mech- 
anism in Bi Q is quantitative rather than qualitative. 

The model - To study the orbital effects of a parallel 
magnetic field on the fluctuations of the SC, we consider 
the microscopic action for a quasi-2D SC film, obtained 
from the BCS Hamiltonian by a Hubbard-Stratanovich 
transformation followed by an expansion around the sad- 
For T < u,Dq 2 < A , this yields 
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FIG. 1: Circles indicate the measured B% in units of = 
(1.6) 2 §q&TcVq I ("fd) Rq I Rq versus the normalized resistance 
per square, R^/Rq, from Ref. @. Here T c is taken from 
an independent measurement. The solid and dashed lines 
mark a fit of Eq. |TJ) to the data, for A(B C ) and A(B = 
0), respectively, using Kq and vq as fitting parameters. Kq 
determines the maximal resistance films that exhibit a SC 
phase. The fitting values are Kq — 0.37, and vq ~ 4.8 x 
10 33 erg _1 cm -3 . Estimates for Kq based on the proliferation 
of vortex loops give 0.46 < Kq < 0.75, see Fig. [2] vq inferred 
from typical carrier densities in amorphous InO is vq w 1.4 X 
10 33 erg~ 1 cm~ 1 . Inset: The measured B c /T c , versus Rq/Rq. 
A transition driven by spin-exchange effects would imply that 
the ratio B c /T c be independent of Ra- Here, we find a large 
scattering of the data. Moreover, in the presence of spin-orbit 
scattering, increasing the disorder (large Ra) would enhance 
spin-orbit scattering, resulting in a larger B c , whereas the 
data of Ref. H shows the opposite tendency. 



where d is the film thickness, v p = -J (3ir/2)DA the 
amplitude velocity, = v/ 7f D A.$(2dV c VQ + 1) the phase 
velocity, £q = itD/8Aq, V c « 2we 2 d is the Fourier trans- 
form of the short range Coulomb interaction due to ex- 
ternal screening, the density of states, D the electronic 
diffusive constant, and the SC order parameter is param- 
eterized as A = Aope 1 ^, with Ao, the mean field solution. 

A parallel magnetic field alters the mean field solu- 
tion, Ao — ► Aq(B) [l2| . and couples to the gradients of 
the phase V</> — -> V</> — —A. We ignore dynamic fluctu- 
ations of the electromagnetic field. This corresponds to 
assuming an infinite penetration depth of the magnetic 
field, and valid in the limit of very thin films. In the 
London gauge the uniform action is 

S = v dA 2 {B) fd 2 rdr {p 2 /2 [hip 2 - l] + p 2 B 2 /^ 2 }(3) 



where H 2 = "1^ ■ Taking p = ap', and choosing a 
to eliminate the third term in Eq. ([3]), we obtain Eq. (|2|) 
with A Q {B) -> A(B) = A (B)atv{-B 2 /H 2 }, and p -> 
p' . Minimizing the uniform action with respect to A = 
A (B)p', yields A MF = A(B) 0. 

Motivated by existing theoretical models that stress 
the role of p hase fluctuations as the cause of the INS 
behavior [laEII, 

we concentrate on the phase action in 
Eq. ©. Taking r -> r/(y/DA (B)), and rescaling the 
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imaginary time z = tAq(B), Eq. ([2]) becomes 



(4) 



where N± = ppd is the number of transverse channels, 



-Kv dD _ 
Ka = — ^ — e 
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_ R Q -2B 2 /H 2 

2R D 



(5) 



is the stiffness coefficient. The stiffness coefficient de- 
termines the action of twisting the phase of the order 
parameter. When Kq ^> 1, the SC phase is rigid. When 
A'o <C f the phase is strongly fluctuating. Hence, there 
exists a critical value Kq that marks the onset of strong 
phase fluctuations, whose exact numerical value depends 
on the details of the transition. From Eq. |J5J), this implies 
a critical parallel magnetic field, given by Eq. {T]). 

Estimate of Kq - While the critical field in Eq. ([1]) de- 
pends on the details of the transition only through Kq, 
this numerical factor becomes increasingly important in 
the limit i?n Rq- Eq. (JTJ shows that Kq determines 
the limiting value of R\j/ Rq for which B c — > 0, and the 
SC phase disappears. To estimate the value of Kq we 
study the microscopic phase action Eq. @. This quan- 
tum action is in the universality class of the anisotropic 
3D XY model. The system described by the 3D XY 
model undergoes a transition between an ordered phase 
(SC) and a strongly fluctuating phase (INS). Different 
mechanisms can drive the system into a strongly fluctu- 
ating phase, including strong gaussian fluctuations and 
the proliferation of topological excitations. Estimates 
based on the Lindeman criterion give Kq up to a numeri- 
cal factor which is usually determined experimentally. A 
more accurate estimate can be done based on the prolif- 
eration of vortex loops, using a 3D generalization of the 
Kosterlitz Thouless scaling procedure 2(| 21 , 22 , 23[ . We 
note that the critical exponents inferred from the 3D XY 
model are consistent with few experimental results Q . 

Previous works calculated the critical stiffness for a 
phase only 3D XY model, Eq. g]) [HIH. The resulting 
Kq implies that films with Rei/Rq > 0.64 should not ex- 
hibit a transition, since B c = 0. Conversely, the data of 
Ref.Hshows that a transition persists for Ru/Rq < 1.31. 
Moreover, the phase only model is applicable to inhomo- 
geneous systems such as Josephson junction arrays. In a 
homogenous system such as a SC film, however, a vortex 
excitation can only occur once the SC amplitude is lo- 
cally suppressed. We solve the flow Eq. ([6|) and (|7|), with 
corrected initial conditions to account for both the phase 
rotation and the amplitude suppression inside the vor- 
tex core, see the discussion following Eq. ([7]) . The critical 
bare physical parameters are the initial conditions that 
flow to the critical point of Eq. ||7J). The microscopic ac- 
tion [Eq. {2])] allows to express these critical parameters 
in terms of measurable quantities, such as Ru and B. 

The large number of perpendicular channels, N± 1, 
in quasi-2D films, generates a strong anisotropy between 



the spatial dimension and the imaginary time dimension, 
Eq. (|4]) . This strong anisotropy in the stiffness coefficient 
introduces crossover scale a' = £oN±. The scaling of the 
phase stiffness and vortex loop fugacity in the 3D XY 
model can be obtained in two regimes. At small distances 
a < a', dominant excitations are found to be rectangular 
loops, cutting single planes (left inset of Fig. [2l 12411 . The 
scaling equations at these scales are quasi-2D [23j 

dK? D /dl = -4n 3 K? D y? D 

dyf D /dl « [4-2^(l + C/(0)iVlV/2)]y^(6) 

where I = In a is the running scale, and Kf D and yf D 
are the quasi-2D renormalized stiffness and fugacity, re- 
spectively. Here 17(0) = E 9 , w 47r /(? 2 + Uj2 / N D is the 
phase propagator, and the sum is cutoff at the effective 
core size that accounts for the crinkling of the vortex 
loops [23( | . At larger distances a > a' the system is no 
longer sensitive to the anisotropy, and dominant excita- 
tions are multiplane vortex loops (right inset of Fig. [2]). 
In this regime, the renormalized Kf, D and yf, D of Eq. (|6|) 
at I' — In [a! /a) arc used as initial conditions for the 
isotropic scaling equations for multiplane loops Q 



dh'i/dl 
dyi/dl 



= K 
= [6 



- 4ir 3 /3 yi K? 
^Kil-xlnK)] yi . 



(7) 



Here x — 0.6 is the self avoiding random walk exponent. 
It accounts for partial cancelation of the Biot-Savart-likc 
interaction in complicated loop geometries (22. 25|. 

We calculate the initial conditions of Eq. ([6]) for a 
homogenous system, with both phase and amplitude 
fluctuations. The bare value of the stiffness coefficient 
is Kl D = K /[l + (2A0)- 1 ] 0, with A'o given by 
Eq. §5$. The bare fugacity of a vortex loop is yQ D = 
cxp {— Sj — S c }, where Sj and S c are the a self rotational 
and core actions, respectively. Previous works calculated 



Sj for a phase only model 



The self rotational ac- 



tion of the smallest rectangular loop was found to be 
Sj = 7T 2 A' /[1 + (2A ) -1 ]. Here we add an estimate for 
S c , calculated from Eq. using a variational method. 
The anisotropy between the spatial dimension and the 
imaginary time dimension introduces two possible exci- 
tations: a vortex, manifested as a rotation of the phase in 
the x — y plane and a phase slip which is the correspond- 
ing rotation in the i-rory-r planes. A vortex loop in 
the 2 + 1 dimensional world is a complicated combination 
of vortex and phase slips segments. We first estimate 
the core action of a vortex/phase slip segment of unit 
volume, using the core size r$ as a variation parameter 
following Ref . [ljj The resulting core action segments are 
Sr = K [y^j6 + 7r/(2AjJ] and 5 C V = K (1 + 3tt) /4. 
In quasi-2D films with N± 3> 1, rectangular loops domi- 
nate [24( . The smallest rectangular loop has £0 unit sides 
of circulation segments in the r direction and 2£o sides of 
circulation pointing in the x — y plane. As a result, the 
bare core action of a rectangular loop is S c = 2S^ + 4S^ S . 
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The bare critical stiffness, Kq(N±), identified as the 
initial condition derived for the phase and amplitude 
model [Eq. ©] that flow to the critical point of Eq. (J7J), 
is plotted in black in Fig. [5] The same quantity cal- 
culated for a phase only model is plotted in gray. The 
corrected initial conditions give 0.46 < Kq(N±) < 0.75, 
for 10 < Nj_ < 10 4 . This implies that for the physi- 
cally realized values of N± a transition could occur up to 
Ro/Rq < 1-08, compared to R\j/Rq < 0.64 found for 
the phase only model. Hence, the estimates for Kq(N±) 
based on the phase and amplitude model are more con- 
sistent with the data of Ref. [H, that show a transition in 
samples with Ru/Rq < 1.31. Our estimates of iS c ne- 
glected the possible crinkling of the vortex loops. These 
would result in a larger core action, and therefore a larger 
range of R\j/Rq that exhibit a transition. 
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FIG. 2: The critical stiffness, Kg, versus the anisotropy pa- 
rameter, NJ 2 = ( p p rf )2 ) calculated using Eqs. (|6} and iff}, 
with the corrected initial conditions including both ampli- 
tude and phase fluctuations (black curve) and the phase only 
initial conditions (gray curve). The critical bare stiffness was 
calculated assuming quasi-2D rectangular loops [U at scales 
a < £N±, and multiplane loops at scales a > £,N±, see inset. 
We find that K§ derived for the phase and amplitude initial 
conditions is more consistent with the experimental data of 
Ref. H, see discussion following Eq.©. 



Summary - We study the previously disregarded or- 
bital effect of a magnetic field applied parallel to a SC 
film. Wc find that the parallel field reduces the phase 
stiffness and leads to strong phase fluctuation at a critical 
magnetic field that depends on the film's sheet resistance, 
Eq. |T]). This prediction does not depend on the details 
of the transition, it allows to experimentally determine 
if spin or orbital effects drive the transition in the paral- 
lel orientation, and it shows that the data of Ref. [8| are 
more consistent with the orbital mechanism, see Fig. [1] 
A quantitative estimate for the transition field, depends 
on the detailed process by which strong phase fluctua- 
tions lead to an INS behavior. In this Letter we consider 
as an example the proliferation of topological excitations 
as a possible mechanism. We map the microscopic action 
of the 2D film onto the 3D XY model. By solving the 
scaling equations derived for the 3D XY model [23J] with 
corrected initial conditions to account for both amplitude 



and phase fluctuations, we get a better agreement with 
the data of Ref. @. Wc note that the phenomenology 
of an INS behavior induced by strong phase fluctuations 
can be generalized to other mechanisms that reduce the 
phase stiffness in a continuous fashion including magnetic 
impurities and changing the thickness, as long as dissi- 
pation in the cores can be neglected. We would like to 
acknowledge useful discussions with E. Altaian, G. Rc- 
facl. Wc thank K. A. Parendo, K. H. Sarwa, A Goldman, 
A. Johansson and D. Shahar for letting us use their data. 
This paper was supported by an ISF and a DIP grant. 
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